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Summary. Numerical simulations of magneto-convection and analysis of solar mag- 
netogram data provide empirical probability density functions (PDFs) for the line-of- 
sight component of the magnetic field. In this paper, we theoretically explore effects 
of several types of PDFs on polarized Zeeman line formation. We also propose com- 
posite PDFs to account for randomness in both field strength and orientation. Such 
PDFs can possibly mimic random fields at the solar surface. 



1 Introduction 



Magneto-convection on the Sun has a size spectrum that spans several or- 
ders of magnitude and develops eddies with sizes much smaller than the 
spatial resolution of current spectro-polarimeters (about 0.2 arcsec or 150 
km at the photospheric level). Thus the Stokes profiles that we observe are 
always averages over space, time, and along the line-of-sight (LOS). This 
suggests that it would be appropriate to characterize the magnetic field 
responsible for spectral line polarization by a p robability density functio n 
(PDF). The case of Gau s sian PDFs was studied byl Polginov fc Pavlovl (Il972h: 



2007.): . Sampoorna et al.l 



Domke fc Pavlovl (|l979l ): lFrisch et al.l (|2005Ll2006a 

(l2008a[l. In this paper we study PDFs determined from observations by 
Stenflo fc Holzreuter 1*2002'. 2003a'lD) and f rom magneto-convection simula- 
tions bylStein fc Nordlund (2006., see also Vogler et al .l2005h .l ISampoorna et al.l 
(|2008bl ) used these PDFs to compute mean solutions in the macro and micro- 



turbulent limits. Macro and micro-turbulence are also referred to as optically 
thick and optically thin limits. In this paper we discuss in detail the mean 
solutions computed for a more general regime of meso-turbulence. 

It is well known that photospheric photon mean free paths (50-100 km 
in the optical) correspond approximately to the sizes of meso-turbulent mag- 
netic e ddies. Calculation of inean Stokes parameters in th is r egime was consid- 
ercd bv'Landi dcgl'Innoccnti (l994VFri sch et all (|2006bl lah: [Carroll fc Staudd 
(|200a.,2005. ): .CarroU fc Kopi (2007,). We use a Kubo- Anderson Process (KAP) 
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with correlation length l/v (where v is number of jumps per unit optical 
depth) and a PDF to characterize the random magnetic field 



th) 

M 



2006bl) . Using KAP in a Milne-Eddington model atmosphere, 



Frisch et al 



Frisch et al 



( 2006b[) deduce explicit expressions for the mean and rnis fluctuations of the 



emergent Stokes parameters. In this paper we use those expressions to com- 
pute mean solutions for magnetic eddies of arbitrary size. 



2 Scalar PDFs 



Recently IStenflo fc Holzreuteil (|2002l l2003allbl ) have found from an analysis 



of high resolution La Palma and MDI solar magnetograms that the PDF for 
the LOS component is nearly independent of the spatial scale and can be 
well represented by a Voigt function. This PDF has a Gaussian core (centered 
around zero field) with a magnetic width Z\b, and Lorcntzian wings with a 
magnetic damping parameter as- The quantity Zi^ is a measure of the rms 
fluctuations of the LOS component. 

A Voigt PDF with Ab = QG and as = 1.5 best fits the empirical PDF 
for LOS field strength derived from La Palma magnetograms. A symmetric 
Voigt function has zero net flux. In real magnetograms magnetic flux is locally 
unbalanced. Non-zero net flux can be generated by a Voigt PDF that has a 
symmetric core, but different as for positive and negative polarities. 



M agneto-convection simulation near the solar surface bv lStein fc Nordlund 
show that the magnetic field is intermittent with a stretched exponen- 
tial distribution. Here we restrict ourselves to asymmetric scalar PDFs as they 
represent the photospheric conditions better. 

2.1 Asymmetric Voigt PDF 

In terms of parameters as and As, symmetric Voigt PDF has the form 

P^iy,as)^^ ^—-^^du. (1) 

In the above equation we have introduced non-dimensional parameters 

V^B/Ab, u^Bi/Ab, jb^Ab/Bj,, (2) 

where 1/Bd = ge / (AnmcAi'Y)) , in standard notation. Here B is the random 
magnetic field component along a given direction, and Ai'd is the frequency 
Doppler width. Thus jb represents rms fluctuations Ab converted to Zeeman 
shift in Doppler width units. 

Asymmetric Voigt PDFs can be constructed by choosing different values 
of as for different parts of the PDF while keeping the Gaus sian core symmet 



rical. Figure [T^ shows examples of asymmetric PDFs (see ISampoorna et al 
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Fig. 1. (a) Asymmetric Voigt and (b) Stretched expo nential PDFs. The y-scaXe i s 
related to the -B-scale through Ab (= 6 G according to lStenflo fc Holzreutejbood ). 



2008bl for details), which more or less resemble the PDF for the La Palma 
magnetogram shown in Fig. 2 of lStenflo fc Holzreuter (l2002h . The mean mag- 
netic field ?/o is the average of y over Pv(y,aB)- 

Mean emergent residual Stokes parameters {rj^q^uy) are computed using 
asymmetric scalar PDFs for magnetic eddies of arbitrary size. The random 
field has a fixed orientation with respect to the LOS defined by polar angles 
(Tj x) = (60°, 30°). Emergent residual Stokes parameters are defined by r/ — 
[Ic - I]/Ci and rx = -X/Ci, with the symbol X denoting Q, U or V. The 
constant Ci is the slope of continuum source function S{tc) — Cq + CiTc, 
with Tc the continuum optical depth. The continuum intensity at the surface 
is /c = Co + Ci. This model has damping parameter a = and line strength 
parameter (ratio of line to continuum absorption coefficient) /3 = fco/fcc = 10. 
We also assume that the spectral line has a wavelength around 5000 A, a 
Lande factor of 2 and a Doppler width of 1.5kms~^. For this typical line, 
Bd = 1.07 X 10'^ G. Hence 7s ~ 0.0056 for an rms magnetic field fiuctuation 
Ab = QG. As a, result the (r/) profiles in all figures in this paper remain 
insensitive to the PDF parameters. 

Figured] shows {ri^q^uy) computed using the 3 PDFs in Fig. [T^ for v — b. 
Note that (ry) is zero if one uses a symmetric PDF with zero mean field. The 
{rq^u) profiles show very small sensitivity to PDF asymmetry, while the (r/) 
profiles indeed are insensitive. For all three PDFs, (ry) peaks around x « 1.5. 
The peak amplitudes increase with the mean field. 



2.2 Asymmetric stretched exponential PDF 

A symmetric stretched exponential may be written in functional form as 

P,e(2/)d2/ = Ce-l^l'd2/. (3) 
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Fig. 2. {ri^Q,u.v) computed with asymmetric Voigt PDFs in Fig. [T^ for = 5. 
Reference to line types is tiie same as in Fig. [1^. 



The quantity k takes values between and 1 and is referred to as the stretching 
parameter. C is the normalization constant. 

We construct asymmetric stretched exponential PDFs with non-zero mean 
field by choosing different k valu es for positive and negat ive polarities. Fig- 
ure [TJa shows three examples (see Sampoorna et al. 2008bl for details). 

Mean residual Stokes profiles computed with these three PDFs are shown 
in Fig. [3l For i/ = 5, we observe a clear increase in the peak amplitudes of 
{fQ,u,v) when the mean field yo increases. However, the peak positions are 
essentially insensitive to j/o- The (r/) profiles remain insensitive to PDF asym- 
metry. Differences between the solutions for = 5 and the macro-turbulent 
limit appear for {rq^u) when yo = 11.6, due to the extended tail of the PDF for 
positive polarities. For (ry), the differences remain small even for yo = 11.6. 
The relative insensitivity of (ry) to the scale of magne tic field fluctuations has 
been observed for Voigt type PDFs (see lSampoorna et al. 2008b) and also for 
symmetric Gaussian PDFs with non-zero mean field (iFrisch et al.ll2006bh . It is 
due to the fact that the PDF asymmetry is sharply peaked around yo- In the 
limit of a Dirac distribution, there would be no difference between solutions 
for different ly values since the magnetic field would be deterministic. It thus 
appears that a mean value of Stokes V can be calculated with rea sonable con- 
fidence using the micro-turbulent limit, a remark made already in lFrisch et al.l 

diooeH). 
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Fig. 3. {ri^Q^u,v) for asymmetric stretched exponential PDFs. Solid, dotted and 
dashed lines refer to v = 5 (the same line types as in Fig. ^p). Heavy dot-dashed 
line refers to macro-turbulent limit for yo = 11.6. 



3 Angular PDFs 

A large fraction of the solar atmosphere is filled with mixed polarity fields, and 
the inter-granular lanes contain fields that are preferentially directed upward 
or downward. To represent this situation, we consider magnetic fields that 
have a fixed value of the strength B but random orientations. For such a 
rando m field, following angular distribution has been suggested by IStenflol 



(1983): 



PMPB)^^^^\f^B\^, -1<MS<+1. (4) 

Here, /is = cos 9b, with 6b the field orientation with respect to the atmo- 
spheric normal. The power law index p can take any value; the p = case 
corresponds to an isotropic distribution. As p increases the distribution be- 
comes more and more peaked in the vertical direction. 

Mean profiles {ti^q) are calculated for z/ = 5, magnetic field strength 
B / Bd =0.1, line strength /3 = 10, damping parameter a = 0, and the power 
law index p a free parameter. Figure IH s hows (r/.p) at the limb (u = 0.1). A 



comparison of this figure with Fig. 12 of lSampoorna et al ( 2008b ) computed 



for the micro-turbulent limit shows that there is very little difference between 
the two. This is because the absolute value of the magnetic field along the LOS 
is bounded by the condition B/B^, =0.1. Since the field is weak, (rg) ^ (rj). 
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Fig. 4. {vi^q) at fj, — 0.1 (limb observation) computed for u = 5 and angular power 
law PDF (Eq. (U)). Line types: p = (solid), 5 (dotted), 10 (dashed), 100 (dot- 
dashed), 500 (dash-triple-dotted), and 1000 (long dashed). In this case {ru,v} ~ 0. 



Further, (rg) is zero for p — 0. As p gets larger, (rg) first increases with p 
and then saturates for p ~ 100 (see Fig. 



4 Vector magnetic field distributions 

To describe a random vector magnetic field we need a PDF that simultane- 
ously accounts for strength and angle fluctuations. From physical considera- 
tions one may argue that angular variation should be strongly field-strength 
dependent (see de Wiin et al. 2009f) . For the strongest fields the distribution 



should be peaked around the vertical direction, as strong fields tend to have 
intermittent fluxtube morphology. The weakest fields on the other hand would 
be passively moved and bent by turbulent fluid motions and so get so tan- 
gled up that their distribution would be nearly isotropic. The transition from 
isotropic to peaked distributions would probably be gradual (possibly around 
50 G). 

Based on this scenario we propose PDFs here which are combinations of 
angle and field strength distributions. For the angular part we use the power 
law distribution introduced in Sect. [H while for the field strength part we 
consider either a Voigt function (see Sect. 12. 1| ) or a stretched exponential (see 
Sect. I2.2p . The functional form of such a composite PDF is 



, N , , , , /Pv(y,as)^^d/isd0sdy, 



27r \C e 1^1 fjFg dfiB d'pB dy. 

Here y varies in the range [— ?/max, -t-J/max], Ms in the range [0,1] and (f>B in 
the range [0,27r]. Py{y,aB) is given in Eq. ([T]). Asymmetric composite PDFs 
can be constructed by choosing different as (or k) values for positive and 
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Fig. 5. {ri^Q,u.v) for = 5 and = 0.1. Composite PDF containing asymmetric 
Voigt function with j/o = 4.5 (dashed line of Fig.[T^) is used. Line types are : yt — oo 
(solid), yt = 50 (dotted), yt — 10 (dashed), and yt ~ 5 (dash-dotted). 



negative polarities of the Voigt (or stretched exponential). The angle and 
strength distributions are coupled by letting the power law index p depend on 
y. We have chosenp = \y\/yt with yt — Bt/As, where Bt marks the transition 
between isotropic and peaked. Note that yt ~ oo corresponds to fully isotropic 
distribution for all field strengths. 

We calculate mean residual Stokes parameters for the composite PDF 
containing an asymmetric Voigt with mean field yo = 4.5 (dashed curve in 
Fig. UK)- The model parameters are (a,/3, 75) = (0,10,0.0056) and /x = 0.1. 
Figure [5] shows the solutions for z/ = 5. The different curve types correspond 
to different yt values, (r/) is insensitive to the PDF asymmetry and to the 
variation of yt due to the very weak value of 75. When yt 00, the PDF 
becomes fully isotropic, and hence (rg) 0. As yt decreases, the PDF be- 
comes more and more anisotropic and hence (rg) as well as (ry) increase in 
magnitude. In this case {rjj) is generated through the magneto-optical effects. 
Therefore (ru) is very small, with a behavior similar to (tq). 



5 Conclusions 



We have presented mean Stokes profiles formed in media having spatially unre- 
solved magnetic structures with sizes that are comparable to photon mean free 
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paths, using PDFs that describe fluctuations of the ambient field. A Gaussian 



PDF with isotropic or anisotropic fluctuations was considered in iFrisch et al 



lere, we experimented with other types of PDFs with restriction to 



asymmetric PDFs which can generate non-zero net flux as diagnosed by the 
shape of (ry) profile. We consider very weak fluctuations of the magnetic field 
(73 = 0.0056). Thus (r/) profiles remain insensitive to the shape of PDF. In 
contrast, other mean Stokes profiles are quite sensitive to the choice of PDF. 

For a complete description of the random vector magnetic field we need 
PDFs which describe both the angular and strength fiuctuations. We con- 
structed such empirical PDFs by combining a power law (for angular distri- 
bution) with a Voigt function or a stretched exponential (for field strength). 
At the solar surface, weak fields are observed to be nearly isotropic and strong 
fields more vertical. Composite PDFs can simulate such a situation. 
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